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Abstract. Periodic orbits in chaotic systems form clusters, whose elements tra- 
verse approximately the same points of the phase space. The distribution of cluster 
sizes depends on the length n of orbits and the parameter p which controls close- 
ness of orbits actions. We show that counting of cluster sizes in the baker's map 
can be turned into a spectral problem for an ensemble of truncated unitary ma- 
trices. Based on the conjecture of the universality for the eigenvalues distribution 
at the spectral edge of these ensembles, we obtain asymptotics of the second mo- 
ment of cluster distribution in a regime where both n and p tend to infinity. This 
result allows us to estimate the average cluster size as a function of the number 
of encounters in periodic orbits. 



1. Introduction 

1.1. Motivation. Periodic orbits have been at core of much of the mathematical 
work on the theory of the classical and quantum dynamical systems. One of the 
primary motivations to consider periodic orbits in Quantum Chaos stems from their 
role as a bridge between two theories. In the semiclassical limit spectral density 
of chaotic Hamiltonian systems can be expressed by the Gutzwiller trace formula 
through the sum over periodic orbits 7, where each term is determined by the stabil- 
ity and the action of 7 [I] . Straightforwardly, this implies that spectral correlations 
between energy levels of quantum chaotic systems must be related to correlations 
between actions of periodic orbits. This remarkable connection has been fruitfully 
employed in a last decade to deduce universal features of quantum spectrum in 
chaotic systems described by Random Matrix Theory [2J [3] . 

On a heuristic level, the action correlation mechanism can be understood in the 
following way. A generic long periodic orbit in a chaotic system has a number of 
encounters where it closely approaches itself in the phase space. It can be shown 
that for any such orbit there exist partner orbits following approximately the same 
path but in a different time-order. The switches in the direction of motion occur at 
the encounters, see fig. QJi. These periodic orbits form families, or clusters as we refer 
to them hereafter. All members of the cluster traverse through approximately the 
same points of the configuration space and, therefore, have close actions. Moreover, 
if the time reversal symmetry is broken the orbits from the same cluster must be 
close to each other also in the phase space of the system. 

The clustering phenomenon can be rigorously described using symbolic dynamics. 
Assuming that dynamical system allows a finite Markov partition any periodic orbit 
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can be represented as a cyclic sequence of symbols from certain alphabet [3]: 

X = X\X 2 ■ ■ ■ X n . 

To simplify the treatment, we will assume in this paper the simplest possible sym- 
bolic dynamics - alphabet of two symbols Xj e {0, 1} with any sequence allowed by 
the trivial grammar rules. Such symbolic dynamics appears for instance in baker's 
map [5j. Two periodic orbits from the same cluster induce two sequences x, x' with 
the property, that each subsequence of p consecutive symbols in x appears the same 
number of times in x', and vice versa. Accordingly, all periodic orbits can be ar- 
ranged into a number of clusters. It is straightforward to see by the above definition 
that two periodic orbits from the same cluster are separated by metric distances 
of the order 2~ p in the phase space [5]. Consequently, the number p controls the 
differences between their actions. 

In chaotic systems the number of periodic orbits grows exponentially with their 
length. Because of this, there are in general many periodic orbits with close actions 
which belong to different clusters. Nevertheless, clustering phenomenon plays very 
important role in Quantum Chaos. Due to their topological nature the action cor- 
relations between orbits within the same cluster are rigid. Namely, periodic orbits 
from the same clusters keep their action differences small also under small pertur- 
bation of the system. It can be expected that the only systematic contribution into 
spectral correlations comes from such orbits. Taking this into account it seems to 
be a natural and important question to ask what is the distribution of cluster sizes 
for given p and n. Note that in a nutshell this is a purely combinatorial problem, as 
we need to count the number of sequences of length n satisfying certain constraint. 

1.2. Formulation of the problem and main results. The cluster distribution 
can be studied through fc-th moments: 

Z k (n) = Y,\C\\ 
c 

where the sum runs over all clusters and \C\ is the number of periodic orbits in 
the cluster C. It worth noting that the first momentum Z\{n) = \C\ provides 
just the total number of periodic orbits of the period n, whose asymptotics for the 
baker's map is given by 2 n /n. In this paper the consideration will be restricted to 
the second moment Z 2 (n). The ratio (\C\) := Z 2 (n)/Zi(n) then can be used as a 
measure of the average cluster size. 

As can be easily understood the asymptotic behavior of Z 2 (n) for long orbits 
crucially depends on the interplay between parameters p and n. In our previous 
work [B] and related studies [TJ. E] the long trajectory limit n — > oo, with p being fixed 
was considered. Note however, that in a proper semiclassical limit the parameters 
n,p are related to each other. Indeed, the spectral correlations on the scale of 
mean level spacing A are determined by periodic orbits with the periods of the 
order of the Heisenberg time HA -1 . For two-dimensional systems this time scale is 
proportional to the inverse Planck's constant fir 1 . On the other hand, the action 
differences between correlating orbits must be of the order of Planck's constant. 
Since action differences and periods of orbits are controlled by p and n respectively, 
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in a semiclassical regime both parameters must tend to infinity. The main focus 
of the present paper is on the limit: n = tyN, where t is a fixed parameter and 
N = 2 P_1 — > oo. Speaking informally, this regime can be interpreted as the limit 
where the average number of encounters in periodic orbits is fixed. This is so, 
since for long periodic orbits with n ^ N the number of encounters is proportional 
to n 2 /N. Note that this regime is essentially different from the aforementioned 
semiclassical limit needed for the evaluation of spectral correlation on the scale of 
mean level spacing. In the last case N ~ n and the number of encounters in periodic 
orbits tends to infinity. 

The central idea of our approach is to represent Z 2 (n) through the spectral form 
factor of an ensemble of truncated 2N x 2N unitary matrices. As we show, the 
relevant spectral information necessary for evaluation of Z 2 {n) depends on the re- 
lationship between n and N. In the semiclassical regime N ~ n the asymptotic 
form of ^(n) is determined solely by the distribution of the largest eigenvalue in 
the ensemble of random matrices. On the contrary, in the regime yN ~ n the 
relevant information arrives from the edge of the spectrum, where both eigenvalues 
distribution and their correlations are of importance. Remarkably, the density and 
correlations of eigenvalues at the 1/y/N vicinity of the spectral edge exhibit uni- 
versal properties. This observation allows us to obtain asymptotics of Z 2 {n) and 
estimate "average" cluster size (\C\) analytically. 

The paper is structured as follows: In Sec. 2 we relate function Z 2 {n) to the 
spectral form factor of an ensemble of truncated 2N x 2N unitary matrices. We 
then analyze the spectrum of these matrices in Sec. 3. We demonstrate numerically 
that at the edge of the spectrum a universality of eigenvalue distribution holds. This 
universal behavior is then analyzed in Sec. 4 using circular ensemble of truncated 
unitary (CUE) matrices with the standard invariant measure. For these ensembles 
we analytically derive asymptotics of the spectral form factor in the regime yN ~ n. 
In Sec. 5 we exploit these results in order to obtain asymptotics of i^(n) an d consider 
implications for the average size of clusters. Discussion of different regimes along 
with open questions for further research are put in Sec. 6. 

2. Cluster counting as spectral problem 

As a first step, we show below that the problem of finding cluster sizes is equivalent 
to the one of counting closed paths of the same length on a special class of graphs 
G p . Specifically, for the symbolic dynamics of backer's map G p happen to be famous 
de Bruijn graphs. 

2.1. Clusters of periodic orbits on graphs. Let X s be the set of all possible 
sequences of zeroes and ones having the length s. The graph G p is constructed in 
the following way. First, with each sequence [aia 2 ■ ■ ■ a p -i] G A p _i we associate a 
vertex of G p . Second, a sequence a = [a\a 2 . . . a p ], a e X p defines the directed edge 
e a of G p which connects vertex [ ] to the vertex [02^3 . . . aJ . Thus the 

total number of edges is 2 P and each vertex has two incoming and two outcoming 
edges, see fig. [Tb. It is worth mentioning that the resulting graphs G p are actually 
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well known in the literature. These are so-called de Bruijn graphs, first introduced 
in [H]. 




FIGURE 1. On the left is shown a caricature of two periodic orbits 
belonging to the same cluster. Apart from the encounters (red rectangular 
boxes) two orbits closely follow each other in the phase space. At the 
encounters orbits switch the directions, as shown in the inset. On the right 
is the graph G p for p = 4. Each sequence [01020304], Oj 6 {0, 1} encodes 
edges of the graph. 

It is straightforward to see that any closed path j x on the graph G p passing 
through n edges can be uniquely represented by a cyclic sequence x = X\X 2 ■ ■ ■ x n . By 
such identification each edge of G p traversed by j x corresponds to a certain segment 
[ximod nii+imod n ■ ■ ■ ^p-i+imod n] of the sequence x. Let n(x) = {n a , a G X p } denote 
a set of integers, such that n a is the number of times the path 7^ passes through 
the edge a. Correspondingly |n(a;)| = ^2,n a = n is the length of the path. Then 
two sequences x,y G X n , belong to the same cluster if and only if n(x) = n(y). 
Equivalently, x, y are in the same cluster if j x and 7^ traverse edges of G p the same 
number of times (which can be also zero). As a result, each cluster of periodic orbits 
can be unambiguously labeled by the vector of integers n = {n a , a G X p }. Keeping 
this in mind we introduce notation C n for the cluster of periodic orbits corresponding 
to vector n and aim to find asymptotics of: 

(2.1) Z 2 (n) = l C -| 2 - 

\n\=n 

Remarks 2.1. By attaching incommensurate "lengths" to edges of the graph G p one 
can turn it into a metric graph. Then each C n can be seen as a cluster of closed paths 
with equal lengths. Accordingly, counting of cluster sizes is equivalent to counting 
of degeneracies in the length spectrum of the corresponding metric graph. The last 
problem has been studied in [TUl El Ell EE2] for some classes of metric graphs. 

Rather then consider Z 2 (n) directly it turns out to be convenient to introduce a 
slightly different function. Recall that each periodic orbit is uniquely identified by 
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a sequence of symbols up to a cyclic shift. This means all n sequences: 

e X n , i = 1 



n 



\pC(i mod n)-E(i+l mod «.)••• •^(i+n—l mod n)., 

correspond to one and the same periodic orbit 7 X . For prime periodic orbits, which 
have no repetitions all these sequences are different. On the other hand for periodic 
orbits with repetitions some of the above sequences coincide and the total number of 
associated sequences is a fraction of n. Let us denote C n the set of all sequences of 
n symbols corresponding to a cluster C n of periodic orbits. In analogy with Z 2 (n), 
Z\{n) we define the following functions: 



(2.2) 



Z 2 (n) = IC-I 2 ' Z 1 {n) = \Cr. 



\n\=n 



\n\=n 



where \C n \ is the number of sequences in C n . For prime periodic orbits we just have 
\C n \ = n\C n \. Since vast majority of periodic orbits are prime[J see e.g., [I], the 
connection between two functions is given by [5J: 



(2.3) 



rr 



rZ^nXl + Oin- 1 )), 



ZAn) 



1 



n 



^(nXl + O^- 1 )). 



2.2. Spectral problem. To find the size \C n \ of n'th cluster we need to count the 
number of closed paths which go through the edges a G X p of G p exactly n a times. 
To this end we introduce the connectivity matrix Q between edges of the graph and 
the auxiliary diagonal matrix A(0): 



(2.4) 
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The dimensions of these matrices are 2N x 2N, with iV := 2 P_1 . 

The usefulness of introduction of matrices Q and A(0) is seen from the relation 
connecting traces of their powers with the sizes of the clusters C n : 



(2.5) Tr(QA(0))" = ^ \C n \ exp (i(n, 0)), (n, 0) = ^ n ^ 



n 



E 



where the first sum runs over all clusters C n . Eq. (12. 5p allows to express the second 
moment Z 2 as 



(2.6) 



Z 2 = J2 l^l 2 =< |Tr(QA(0))"| 2 >4n 



1 In particular, if n is a prime number there are only two non-prime periodic orbits. They 
correspond to sequences of all zeroes or ones, respectively 
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with the average taken over the flat (Lebesgue) measure: 

To simplify notation we denote the product of Q and A(</>) by 

Q(<f>) := l -QK{<t>). 

Since half of the eigenvalues of Q are zeroes [13] . the matrix Q(4>) has N zero 
eigenvalues for a general <fi, with the remaining N distributed between and 1. By 
eq. (I2.6p the function Z2 can be expressed through the form factor of the ensemble 
of matrices Q{4>): 



2 2n 1 N 

(2-7) Z 2 = — K(n,N), K(n, N) = - < \ £ K{<t>)? >«, 

i=l 

where Aj(0), i — 1, . . . N are non-trivial eigenvalues of Q{4>). 

The representation f !2.7p is the key component of our analysis, as it allows to 
relate Z2 to the information on eigenvalue distribution in the ensemble of matrices 
Q{4>). For forthcoming treatment it is convenient to split K(n,N) into diagonal 
and non-diagonal parts: 
(2.8) 

TV N 

K^ d \n,N) = IW)| 2n K^\n,N) = - < (</>) >^> ■ 

i=l ' ' j^j 

The diagonal part is determined solely by the density of the eigenvalues: 

1 - 

i=i 

Since the flat measure is invariant under the shift (pi — > (pi + c, i — 1, . . . 2iV ', p(z) 
depends only on the modulus a; = \z\ 2 . Taking this into account the diagonal part 
of the form factor can be written down as: 

(2.9) K {d \n,N)=7i [ p{x)x n dx. 

Jo 

In contrast, the non-diagonal part depends on the correlations between eigenvalues. 

3. Ensembles of truncated unitary matrices 

It is a simple observation that Q(4>) are spectrally equivalent to the matrices of 
the truncated unitary type. More specifically, we can find a matrix [7 , such that 

(3.1) U Q(4>)Ul = PU(cj>), P = diag {1, . . . , 1, 0, . . . , 0} 
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with t/(0) being a unitary matrix and P is a diagonal matrix containing iV zeroes 
and N ones. Indeed, by using the unitary matrices 



(3.2) 



( 1 




\0 



1 
1 









1 

-1 







Inxn 
-Inxn 



Inxn 
Inxn 



where Inxn is the N x N unit matrix, we get the above representation: 

(3.3) U Q{<t>)Ul = P{7(0), {7(0) = Uik{<t>)Ul 

The spectrum of matrices PU, where U is drawn from the circular unitary en- 
semble (CUE) has been previously studied in [TH] for the projection P of a general 
rank. In particular, if, as in our case, the rank of P is half of its dimension the 
limiting iV — > oo density of eigenvalues in these ensembles is given by: 

(3.4) p Q (x) : = lim p^^x) = —r— — for x < 1/2 , 

AT->oo 7T(1 — X) A 

and po(x) = otherwise, see fig. [2] The circle \z\ — I/a/2, where the limiting density 
abruptly goes to zero is referred as the edge of the spectrum. For a finite N the 
edge has a finite "width" of the order 1/yN. The uniform asymptotics of spectral 
density at the vicinity of the edge has been found in [T7]. For the particular form 
of the projection matrix P it reads: 

(3.5) PtiVE (x) = - 1 (l - leric (y/AN {1/2 - x)) ) + O(N^), 



7T(1 — X) 



erfc(y) = —= 



~ x2 dx. 




FIGURE 2. The spectral density for matrices Q(0) in comparison to the 
spectral density of truncated CUE (solid blue line). The red (dashed) line 
shows limiting eigenvalue distribution (|3.4p . Systematic deviations can be 
seen in the bulk of the spectrum. 

The flat measure for the ensemble of matrices <2(0) obviously differs from the 
invariant measure of CUE. Nevertheless, the numerical simulation demonstrate uni- 
versality of spectrum properties at the vicinity of the edge. Namely, both eigenvalue 
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density and the eigenvalue correlations in ensemble of matrices Q((p) are the same 
as in the ensemble of truncated CUE on the scales l/y/N around the edge. Below 
we summarize these findings in the form of the conjecture: 

Conjecture 3.1. 1) Local density p(x) of eigenvalues Aj(0) is universal on the 
scales 1/y/N around 1/2 i.e., for a fixed s 

(3.6) A„ W := p Q - ^) - p™, g - ^) = o(*ri), 

where p TlV¥ ,{x) is the density of eigenvalues in the ensemble of truncated unitary 
matrices with the invariant measure, whose asymptotics is given by $3. 5\) . 
2)Letn = ty/N and let 

Kit) := 2 n VNK(tVN, N), /C TrUE (t) := 2 n VNK TiVE (tVN, N) 

be rescaled form factors of respective ensembles of truncated unitary matrices^ In 
the limit, where t is fixed and N — > oo both spectral form factors have the same 
asymptotics: 

(3.7) A/C(t) := K{t) - /C TrUE (t) = O(N^). 




FIGURE 3. Convergence of the radial eigenvalue density function of 
the matrix ensemble Q((p) to the density of the corresponding truncated 
unitary ensemble in the vicinity of the spectrum edge. The left plot shows 
Ap(s), see (|3.6p . for two matrix sizes. On the right picture the values of 
log 2 < |Ap| > are plotted for different p. The averaging is taken over the 
region s E (—1,1). The straight dashed (red) line with the slope 1/2 is 
shown here for comparison. 

To demonstrate the validity of this conjecture we have numerically calculated 
spectrum of matrices Q((f>) for thousands of realizations of (p. The averaged over 
cp form factor and density of states are shown on figs. [31H1 for different values of N. 
These plots clearly demonstrate convergence of /C(t) and p(l/2 — s/y/N) to their 
counterparts in the corresponding ensembles of truncated CUE. 

2 The rescaling by factor 2" is equivalent to the multiplication of each eigenvalue Xi((f>) by \/2. 
Note that this shifts the edge of the spectrum of matrices Q(<f>) to 1. 
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FIGURE 4. Rescaled spectral form factor /C TrUE (t) of truncated CUE 
is shown in the figure (a) as function of the parameter t = n/y/~N for 
several parameters p (dashed lines). In the figure (b) is similar plot for the 
rescaled spectral form factor lC(t) of the ensemble of matrices Q{4>)- The 
red solid lines depicts asymptotics f sinh(i 2 /4). On the bottom figures are 
the plots of the difference A/C(t) as functions of t (c) andp (d), respectively. 
The dashed straight line with the slope 1/2 is depicted in figure (d) for 
comparison. 



Remarks 3.2. 1 ) Since the diagonal part of the form factor is determined by density 
p only, the above conjecture actually implies that both K^, K^ nd > asymptotically 
converge to their CUE counterparts. 

2) Recall that in the present work we consider symbolic dynamics arising in the 
baker's map. For more complicated grammar rules the matrix Q((f>) is not necessarily 
spectrally equivalent to a matrix of truncated unitary type. In general this is a sub- 
unitary matrix which can be brought (by some unitary transformation ) to the form 
DU , where D is a diagonal matrix with elements in the interval [0, 1] and U is a 
unitary one. We believe that universality at the edge of the spectrum holds for these 
matrices as well. Numerical simulations show that spectral density and form factor 
asymptotics at the vicinity of the edge can be deduced by assuming that U belongs 
to CUE. The spectral density for such ensembles has been studied analytically in 

3) So far justification of the above conjecture comes mainly from the numerics. 
We believe, however, that analytical methods based on the supersymmetry approach 
[15] can be used to prove the conjecture and leave it for future investigation. 
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4. Uniform asymptotics of the form factor 

According to the above conjecture, in the regime where lengths of periodic orbits 
are of the order y/~N one can use the spectral form factor of truncated CUE in order 
to extract asymptotics of Z 2 {n). The goal of this section is to derive the asymptotics 
of K TrVE (n) in the limit, where the ratio n/y/N is fixed. 

It is instructive to separate ii' Tl . UE (n, N) = K^ iVE {n, N) + K^f UE (n, N) into diagonal 
and non-diagonal parts: 

TV N 

(4.1) K^Un) = ^ < E N 2 " >> K ^ n ) = n< >, 

i=l ' ij^j 

where we dropped label N from the definition of K^l E , K^ E for the sake of notation 
compactness. The exact analytical expressions for both parts can be easily obtained 
from the representation 
(4.2) 

KSUn) = d 2 zR 1 (z)\z\ 2n , Jc££(n) = ±J d 2 z x J d 2 z 2 R 2 {z x , z 2 )(z h z* 2 ) n , 

where R\, R 2 are one-point and two-point spectral correlation functions, respectively. 
Using the results of [TB] for R ly R 2 we obtain for the diagonal part: 



(4-3) K%L(n) = ^fl 



' F(m + n)T(m + N) 



in- 



? T(m)T(m + n + N) 



with r(s) standing here for the gamma function. The non-diagonal part is given by 



(4.4) K^ E {n) = -i £ 



A n T(m + n)T(m + N) 



N ^ T(m)T(m + n + N) 

m=l v ' v ' 



for n < N, 



and K^ E (n) = otherwise. To proceed further we use integral representation of 
gamma functions, 

1 + e iv ) N+m ~ l e-^ N dv 



r(m + n + N) 
T(m + n) N f2n 



T(N)T(m) 2vr 
to convert sums (14.3|4.4|) into double integrals: 
(4.5) 

rXd) /n f 1 f 2K (l-x) N - 1 X n + N \l + ^>\ 2N J . I 112 , . nj 

K^ E (n) = - I I _ - d V dx + nJ^ p (x)x dx, 



o Jo 



(4.6) 

/o Jo 2n{l-{l + e^)x) 



rf{nd) / s I 1 [ 2 * {l-x) N -'x N \l + e^\ 2N {l + e^)-^ f 1 ' 2 
&Trvk\n) = / / — — — — — dipdx — Tr / po{x)x dx, 
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where po(x) is the limiting spectral density (13 .4p . We will derive now the asymptotics 
of both diagonal and non-diagonal parts when n = ty/N, N — y oo and t is fixed. 

4.1. Non-diagonal part. Let us first analyze the non-diagonal part. After change 
of the integration parameters x — > l/(x + 1), e lLp — > z the integral (14.41) can be cast 
into the form: 

(4.7) 

Tf(nd) r,r x N ~\x + ir^n + zr-z- N -\ W 

K-TrVTs(n) = / ax (b r az — hi Po(x)x ax, 

Jo Jc 2n(x-z) J 

where the integration over z variable is performed along the circle C of the unit 
radius centered at the origin of the complex plane. The asymptotics of the first 
integral I in eq. (14.71) can be found using saddle point approximation. Writing dawn 
the above integral as: 

roc r e NS(x,z) 

(4.8) 1=1 dx idz 



2tt(x — z) 



(4.9) S = ± log (x N ~\x + 1)- 2 ^(1 + zf N -z- N -') , 

one finds that the saddle point of S is at x = l,z = 1 + 0(1/ yN). Note that 
the denominator is of the order 1/y/N at the saddle point. We therefore make 
substitution z = 1 + w, x = 1 + y and expand S up to the quadratic order in w, y, 
keeping the denominator intact. This gives: 

i r+oo f p -nw/2+Nw 2 /4 / i \ 

(4.10) « = « S+ 0(^ 



After rescaling yNw — > iv, yNy — > u and deforming the contour of the integration 
we have: 

T'(i\ / 1 \ 1 r+oo r+oo -« 2 /4— tiv/2— v 2 /4 

(4.11) TI = -^L + O (4t ) , I'(t) = — / dvdu. 



The last integral can be easily evaluated by switching to polar coordinates v 
r cos 6, u = r sin#: 



i r+oo p2n 

(4.12) I'(t) = — / e- r2/A - irtsine/2+ie drd8 

27T Jo Jo 



e- r2/4 J 1 (rt/2)dr = -(1 - e" t2/4 ). 

It remains to find asymptotics of the second term 77 in eq. (14. 7p . The main 
contribution here comes from the boundary 1/2 of the integration. As a result, to 
the leading order of iV we have: 



(4.13) 2"// = ^ l + O - =— + 0, v 



12 



Combining / and // we obtain for the non-diagonal part: 
(4.14) JC^l(t) := y/NFKg&in) = - V' 2/4 + O 



4.2. Diagonal part. The asymptotics of the diagonal part can be obtained in a 
similar way. Alternatively, one can use uniform asymptotics (13. 5p for the density of 
states and plug it into eq. (12. 9p . The resulting expression is again the exponential 
function but with the opposite sign: 

(4.15) jd£L(t) := VN2 n K^Un) = V/ 4 + O (^=) . 

Remarkably, the leading asymptotics of the diagonal and non-diagonal parts turn 
out to be connected by a simple relationship: 

(4.16) iK^(it) ~ K%SL( t y 

5. Estimation of cluster sizes 



The conjecture (13. 1 p implies that for n = ty N the asymptotics of K(n, N) can be 
extracted from eqs. ( 14.151 14. 14j) . Summing up the diagonal and non-diagonal parts 
of the form factor we obtain by eq. ( 12. 7ft : 

(5.1) (^)z 2 =isi„h( t 74) + (J V-J). 

This is the main result of the present paper. To better clarify its meaning it is 
informative to consider an average size of clusters rather than Z^in) itself. Since 
Z\(n) = 2 n , by eq. (12. 3p the average size of clusters in the regime n ~ \/N can be 
estimated as 

(5.2) < \C\ >= ^ =JC f, + O(N^) = 1 sinh(t 2 /4) + O(N^), 

see fig. Let us show now that the parameter t 2 = n 2 /N actually determines the 
number of encounters v in periodic orbits of the length n. To count the number of 
encounters we will follow the original approach of [2]. 

A periodic sequence x = x\%i . . . x n is said to have the encounter a of the length 
s if the subsequence a = [a\ . . . a s ] appears (at least) twice in x i.e.,: 

(5.3) X = X\ . . . Xk 0>i . . . 3?fc_|_g_|_i . . . X m Qi ■ ■ ■ a s ^ X m -j- 5 _j-i . . . x n , 

encounter encounter 

where x k ^ x m ,x k+s+1 ^ x m+s+1 . 

Let [x], [y] be two (non-cyclic) sequences of the length n. We define the function 5^, 
such that 5^([x], [y]) — 1 if first p symbols of [x], [y] coincide and [y]) = 0, 

otherwise. The number of encounters in a periodic sequence x can be determined 
through the function: 

n n 

(5.4) A4(p) = ^^5 p (^[x],^[x]), 

i=l j 
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FIGURE 5. Comparison of RMT result (blue) (|5.2|) with the numerically 
calculated log < |C| > (red) as function of n for p = 9. 

where a is the shift map whose action on sequences reads as: 

a ■ [x x x 2 ...x n } = [x n x 1 . . . x n _i]. 

Whenever the length s of an encounter is larger than p the function M x {p) counts 
it with the multiplicity 2(s — p + 1). One can, however, get read of this redundancy 
by considering the difference: 

(5.5) v{x) = (M x (p) - M x (p + l))/2, 

where each encounter is counted only once. Note that (15. 5p . actually, counts cor- 
rectly only encounters appearing twice in periodic sequences. This, however, does 
not affect the final result, since encounters of larger multiplicities are rare in the 
considered regime. Indeed, for n ^ N the probability of getting encounter fc-times 
scales as n 2 /N k . 

We can now estimate M x {p) applying the uniformity principle. According to it 
long periodic orbits behave in the same way as trajectories with "almost all" initial 
conditions. Applying the ergodicity theorem to the N x {p)/n 2 and substituting the 
"time" average in (15.41) by the average over the phase space we obtain J\f x (p) ~ n 2 /N 
to the leading order in n. By eq. (15. 5 p this gives for the average number of encounters: 

(5.6) v ~ n 2 /AN, n -»■ oo, n<^N. 

Taking this into account we can finally express (15. 7p through the average number of 
encounters: 

(5.7) < \C\ >= isinh u + 0(N~^). 

6. Discussion 

It is well known that spectral correlations in quantum chaotic systems essentially 
depend on considered time scale. The spectral form factor of chaotic systems with 
broken time reversal symmetry has discontinuity in its first derivative at the Heisen- 
berg time. From the semiclassical point of view such a jump can be attributed 
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to a different behavior of periodic orbits for n <ti N and n 3> N. It is therefore 
interesting to see how this transition affects the average size < \C\ > of clusters. 

Eq. (15.71) demonstrates that for orbits with periods of the order of square root 
of the Heisenberg time, the asymptotics of < \C\ > solely depends on the average 
number of encounters u. For n < \/~N, when encounters are rare, we obtain log < 
\C\ >~ 0. In other words, in this regime vast majority of clusters consist of just one 
periodic orbit while the number of clusters is equal to the number of periodic orbits. 
As the average number of encounters v becomes of the order one, an exponential 
growth of cluster sizes starts. For \/N <C n <C N we have by (15. 7p log < \C\ >~ n 2 , 
implying very fast growth of clusters, see fig. [5j 

The above result should be compared with one from [6] for the regime of very 
long trajectories n ^> N - far beyond the Heisenberg time. In that case the whole 
average (12. 7p is dominated by the largest eigenvalue A max (0) at the vicinity of 4> = 
leading to: 

2 n ( N \ N/i 

(6.1) <|C|>= ^(2^j + 

Here the average size of clusters grows with n in the same way as the number of 
periodic orbits (2 n /n) while the number of clusters grows only algebraically with n. 
This is so, since at n > iV almost every point of a generic periodic orbit belongs 
already to some encounter. The further growth of n does not lead to an essential 
increase in the number of encounters. This results in a slower growth of cluster sizes: 
log < \C\ >~ n. 

The transition between the aforementioned regimes is expected at n ~ N. At such 
scales only the diagonal part of the form factor contributes to Z^in). Furthermore, 
the dominant contribution to the integral (I2.9p comes from the interior of the interval 
[1/2,1], where the spectral density is non-universal. Consequently, the results for 
truncated CUE with the invariant measure are not applicable in this case. In general, 
one can assume that for x > 1/2 the density is exponentially decaying function of 
N: 

(6.2) p(x) ~ exp(-JV$(x)), 

where $(x) is a monotonic function of x with $(1) = oo. Substituting this into 
eq. (12. 9p and performing saddle point approximation gives: 

(6.3) < \C\ >~ exp n (~$( Xo ) + log(2a;o)^ , 

where n = tN and x$ is solution of the saddle point equation: 

(6.4) Xo$'(x ) = r. 

Note that at the limit t > 1 (deep beyond the "Heisenberg time") xq ~ 1 and 
< \C\ >~ 2 n in agreement with (16.11) . For a more precise information on the 
transition from the regime (15. 7p to (16. ip one would need to know exact form of the 
function $(x) i.e., the non-universal tail of the spectral density. 
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In conclusion, let us mention few possible extensions of the present paper results. 
First, general symbolic dynamics with finite grammar rules can be treated in a 
similar way. In this case cluster distribution can be obtained from the spectral 
form factor of more general ensembles of sub- unitary matrices investigated in [17] . 
[18]. We expect that the final result (15. 7p stays the same, while connection between 
v and t depends on symbolic dynamics. In particular, for trivial grammar rules 
with / symbols this connection reads as v ~ n 2 {l — 1)/21N [15]. Second, using 
our approach systems with the time reversal symmetry can be treated, as well. 
In this case one should first redefine the notion of cluster. If the time reversal 
symmetry is present, two periodic orbits belong to the same cluster iff they traverse 
approximately the same points of the phase space up to a switch in the momentum 
direction. This additional freedom results in a certain constrain on the elements 
of matrix A(0). Namely, for any sequence a = [ai,...a p ] and its time reversal 
counterpart a* = [a p , ...ai] we have <p a = <p a *. In its own turn this symmetry 
implies that effectively the matrix ^QA(<f>) belongs to the ensemble of truncated 
orthogonal, rather than general unitary matrices. Note that such ensembles with 
the Haar invariant measure were previously studied in [TjJ]. We believe that these 
results can be exploited in order to obtain information on the cluster distribution 
for periodic orbits in the dynamical system with time-reversal invariance. 
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